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The mechanical response of a wet granular layer to imposed shear is studied experimentally at low 
applied normal stress. The granular material is immersed in water and the shear is applied by sliding 
a plate resting on the upper surface of the layer. We monitor simultaneously the horizontal and 
the vertical displacements of the plate to submicron accuracy with millisecond time resolution. The 
relations between the plate displacement, the dilation of the layer and the measured frictional force 
are analyzed in detail. When slip begins, the dilation increases exponentially over a slip distance 
comparable to the particle radius. We find that the total dilation and the steady state frictional 
force do not depend on the driving velocity, but do depend linearly on the applied normal stress. 
\Q • The frictional force also depends linearly on the dilation rate (rather than the dilation itself), and 

reaches a maximum value during the transient acceleration. We find that the layer can temporarily 
sustain a shear stress that is in excess of the critical value that will eventually lead to slip. We 
describe an empirical model that describes much of what we observe. This model differs in some 
respects from those used previously at stresses 10 6 times larger. 

(N . PACS: 83.70.Fn, 47.55.Kf, 62.40.+i, 81.40.Pq 
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I. INTRODUCTION 



The response of any material to an applied shear stress is an important mechanical property. We are concerned 
here with granular materials, for which the response to shear stress amounts to friction. We have previously reported 
an extensive study of dry granular friction under low applied normal stress and here we extend that work to 

the case of wet materials. Although our motivation is primarily fundamental, it should be recognized that a study of 
granular friction, including the wet case, is of potential interest in connection with applications such as the processing 
of powders. It may also be of conceptual or heuristic interest in connection with seismic phenomena, though the 
applied stresses in our work are lower by a factor of 10 6 . Do the friction laws obtained at pressures of 100 MPa 
extrapolate to such low stresses, or do completely new effects occur? An excellent review of geophysically inspired 
experiments conducted at high pressures has been given by Marone || . Much of the work described in that review was 
concerned with the determination of friction laws that incorporate memory effects. These state variable friction laws 
are also reviewed in the book by Scholz p| , and a recent review of friction that includes consideration of microscopic 
physics has been given by Persson 0. 

We have several different motivations in considering the case of a granular medium that is immersed in a fluid. 
It is well known that even slight amounts of water adsorbed on a granular medium (or a rough surface) can affect 
friction. In one recent paper, the condensation of small amounts of water near points of contact was shown to produce 
substantial strengthening of a granular material ||. To eliminate this effect, one can fill the pore spaces entirely with 
a fluid, i.e., submerge the material. Once we undertook such an experiment, we found that immersion has additional 
merits : Lubrication reduces the friction so that the transition from stick-slip dynamics to continuous sliding occurs 
at very low velocity. Thus, one can study the dynamics of the transition from rest to steady sliding over a wide range 
of velocities, which turns out to allow delicate tests of friction laws. 

Broad references to the literature on granular friction were given in the introduction to our previous paper [[lj, 
and we do not repeat that material here. Works specifically concerned with friction in wet granular materials appear 
to be rather limited. However, it is important to mention the experimental work of Marone et al. on the frictional 
behavior of simulated fault gouge 0. They studied the response of immersed granular material subjected to shear 
deformation and observed both a noticeable dilation of the material and a consequent velocity strengthening. These 
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results were obtained at large normal stress (50-190 MPa) in sand. The work presented here differs from their work in 
two essential respects: First, the normal stress applied to the layer is at least six orders of magnitude smaller in our 
work (typically 20 Pa). Second, the granular material considered here consists of spherical glass beads with a narrow 
size distribution. We point out in a later section the similiarities and differences in the resulting frictional dynamics. 
A thorough controlled study of the frictional behavior of an immersed granular material under low stress, comparable 
to what we have presented previously for dry materials, has not appeared previously to our knowledge. 

The experimental setup is similar to that described in Refs. The sensitivity achieved in the present paper 

is an advance over previous work. We are able to measure both vertical and horizontal displacements to a precision of 
about 0.1 /im, and are able to do so with excellent time resolution of < 0.1 ms. It is true but surprising that precision 
for displacement measurements much better than the particle size is necessary in order to explore the dynamics fully. 
We focus attention here particularly on the vertical dilation that accompanies slip, since it is clear that dilation is a 
necessary dynamical variable in any granular friction law. 

The present article is organized as follows. Our experimental methods are described in Sec. II and the observations 
are presented in detail in Sec. III. In Sec. IV we compare our results to an empirical model and previous work and 
then we conclude. 



II. EXPERIMENTAL METHODS 
A. Introduction 

A schematic diagram of the experimental setup is shown in Fig. 0. The granular material resides in a rectangular 
tray whose lower surface has been coated permanently with a layer of particles to ensure that no sliding occurs at 
this interface when the layer is sheared. A given amount of granular material is then spread out on the tray and it 
is filled with distilled water. A homogeneous layer of material 3 mm deep is obtained by sliding a straight rod along 
the side edges of the tray. The sliding plate is gently placed on top of the granular layer, additional water is added 
to ensure that it is submerged and covered by about 1 mm of water, and a visual check is made to ensure that no air 
bubble has been trapped underneath the plate. 

Shear stress is imposed on the granular layer by translating the plate over its upper surface. In order to transmit 
the shear stress, the bottom surface of the sliding plate is roughened in one of two ways: (a) If visualization is not 
required, a glass plate is used with a monolayer of granular material glued to its lower surface, (b) When visualization 
is desired, a transparent acrylic plate is used; its lower surface is ruled with parallel grooves 2 mm apart and about 0.1 
mm deep to pin the material. Most of the experiments reported here have been done using (b), but we have checked 
that (a) gives essentially the same results. 

The dimensions of the sliding plate are length L = 8.15 cm, width W = 5.28 cm, and thickness T — 0.88 cm; the 
direction of travel is along the length. The in-plane dimensions of the tray (i.e. 11 x 18.5 cm 2 ) are much larger than 
those of the sliding plate; no boundary effect has been observed. 

The granular material used in most of the experiments reported here consists of spherical glass beads with a mean 
diameter of 103 /im and a standard deviation that is 14 /iin from Jaygo, Inc. We have checked that our conclusions 
remain valid for other samples with mean diameters of about 200 and 500 /im. 

Whether or not a layer of beads was glued to the bottom of the container did not affect the experimental results 
described in this paper. We conclude that the shear zone is smaller than the layer height in all cases considered. The 
small total vertical dilation of the granular layer during shear, described in Sec. Ill, indicates that the shear zone 
height is, at most, a few beads. The size of the shear zone can become comparable to the height of the granular 
layer for the largest beads (for which the layer height is approximately six bead diameters) at low pulling velocities. 
However, such a case is not included in the experimental results. 

Visual inspection of the top layer of beads underneath the transparent acrylic plate shows that the beads are 
arranged randomly. The arrangement remains random after repeated sliding of the plate. The organization of beads 
in deeper layers cannot be determined visually. 

The volume fraction of glass beads is 0.63 ± 0.02 for both wet and dry cases. We find that shaking and pressing 
of the material does not increase the volume fraction measurably. The measured value is close to the numerically 
obained estimate of 0.63 for random close packing of uniform spheres p0[ ; a slightly larger value would be expected 
for random close packing if there is some variation in bead size. 
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B. Driving system and frictional force measurement 

In order to control the stiffness of the driving system, the plate is pushed with a blue tempered steel leaf spring 
(constant k — 189.5 N/m) connected to a translation stage. The stage is driven at constant velocity V (between 
0.1 fxxsx/s and 1 mm/s) by a computer-controlled stepping-motor which turns a precision micrometer. The coupling 
between the spring and the plate is accomplished through a rounded tip which is glued to one end of the plate. This 
makes it possible for the plate to move both horizontally and vertically as the thickness of the granular layer varies 
during the motion. 

The elastic coupling through the spring allows relative motion of the plate with respect to the translator. The 
horizontal relative motion is monitored by measuring the displacement d(t) of the spring from its rest position at the 
coupling point with an inductive sensor, model EMD1050 from Electro Corporation, as in Refs. A typical trace 

of d as a function of time is shown in Fig. ^(a) . 

The differential equation governing the position x(t) of the sliding plate in the laboratory frame reads 

= k(Vt - x) - F = kd - F, (1) 

dt z 

where M is the mass of the sliding plate, k is the spring constant, and V is the pulling velocity. Here F denotes the 
frictional force, which depends on time. For the parameter range used in the experiment, the inertia of the sliding 
plate can be neglected and the frictional force is then F — kd. We estimate that the displacement d(t) is known 
within 0.1 fim so that the frictional force F is measured within 2 x 10~ 5 N. This is typically about 0.01% of the 
normal stress. 



C. Vertical displacement measurement 

The vertical displacement h(t) of the plate is monitored with a second inductive sensor as shown in Fig. |||(b). The 
target of the sensor is a blue tempered steel plate connected to the sliding plate with the help of five pillars (diameter 
= 3.4 mm). The target is not perfectly flat and parallel to the sliding plate; the measurements must be corrected 
to exclude this source of error at the micron level of sensitivity. Once this is done, we estimate that the vertical 
position of the plate is known within about 0.2 /xm. We also point out that h(t) is determined only to within an 
additive constant; We choose the constant so that h[t — ► oo) = in the steady-state regime [e.g. the right hand side 
of Fig. |(b)]. 



D. Applied normal stress: effect of immersion 

The plate is immersed in water, so it experiences a buoyant force. The mass to of the displaced volume of water 
(volume of the plate + volume of 1 mm length of the pillars) is to = (12.5 ±0.5) g. In the following, we define the mass 
of the sliding plate M to be its real mass reduced by m. For instance, most of the results given below are obtained 
with a sliding plate whose actual mass is 27 g, so that M = (14.5 ± 0.5) g. The mass of the plate must be small 
enough that the plate does not sink significantly into the granular layer. We estimate the density of the water-glass 
beads mixture to be about 2 g/cm 3 . Then, assuming that the granular layer behaves like a dense fluid, estimating 
the equilibrium depth of the bottom surface of the plate, and assuming that the plate sinks (during translation) when 
this depth equals the plate thickness, we determine that the plate should be expected to sink into the granular bed 
for masses larger than 40 g. Experimentally, we find that the granular layer can sustain somewhat larger loads but 
we limited our study to M < 40 g. 



III. EXPERIMENTAL RESULTS 
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A. Basic observations 



We always prepare the system for measurements in the same way to assure reproducible experimental conditions. 
The granular layer is first installed as previously described, and the plate is then pushed at constant velocity (20 /xm/s) 
until a steady state regime is reached. During this initial "break-in" phase, the plate slightly sinks in the granular 
layer and a small heap forms at the leading edge. This heap has been proven to have a negligible effect on the frictional 
force by turning the plate through 90 degrees. The results obtained by pushing the plate either along its length or 
along its width differ only slightly, despite the factor of 1.5 difference in dimensions. The spring is then pulled back 
until it loses contact with the rounded tip so that no stress is initially applied for the measurements of interest. 

Once the sample has been prepared as just described, the spring is pushed at constant velocity, and both the spring 
horizontal displacement d(t), and the vertical position of the plate h(t) are monitored during the motion as shown in 
Fig. ||. The horizontal position x(t) of the plate in the laboratory frame is deduced from d(t) through the relation 
x{t) — Vt — d(t) (Fig. ||). After a transient regime (which is discussed in Sec. III.C), both the frictional force and 
the vertical position of the plate are found to reach steady asymptotic values, provided that the spring is pushed at 
a velocity larger than about 1 //m/s. We focus next on the behavior in this steady-state regime. 



B. Steady-state regime 

The frictional force Fd in the steady-state dynamical regime is proportional to the mass M of the sliding plate as 
shown in Fig. ^. The frictional force is then written: 

F d = puMg, (2) 

where [id is the frictional coefficient in the steady-state regime, which might depend on the plate velocity. There is 
no evidence of an additional viscous contribution to the frictional force, which could arise due to the surrounding 
water. Indeed, the linear interpolation leads to Fj, ~ for M = 0. [We find Fd = 0.0002 N for M — 0, a value that is 
consistent with zero, given the uncertainty in M .] It is important to know that such a contribution is insignificant, 
since it would be independent of the normal applied stress. The experimental slope in Fig. ^ leads to a friction 
coefficient fid = 0.236 ± 0.004 for the 103 /j,m diameter beads. This value is significantly smaller than that measured 
for the same material when dry. In that case, we observe stick-slip motion for which /id is found to vary between 
about 0.4 and 0.6 §. 

We also analyzed the dependence of the frictional force on the velocity in the steady-state regime that arises after 
the initial transient. Within the experimental resolution, we find that the frictional force does not depend on the 
pulling velocity V over four orders of magnitude, as shown in Fig. |^. 

However, when the velocity is very small (less than about 0.1 /im/s), we find stick-slip motion pd| (Fig. ^|). In 
this case, which corresponds to the leftmost point in Fig. ||, we have used the time average frictional force as an 
estimate. The steady state, which could presumably be reached for such small velocities by increasing the spring 
constant, would probably have a somewhat higher friction coefficient since velocity weakening is a requirement for 
the stick-slip instability. Except for such very low velocities (typically V < 0.1 /im/s), it is safe to conclude that the 
friction coefficient is independent of velocity. 



C. Transient regime 

1. Description 

The sequence of events in the transient regime is as follows. When the spring first comes into contact with the 
rounded tip at A in Fig. g(a), the applied stress increases linearly with time, and the plate remains at rest (Fig. ||) 
until a minimum horizontal applied stress is reached at B. During this initial stress loading over the interval AB, no 
significant change in the vertical position of the plate is observed [Fig. ||(b)]. When the applied stress is large enough, 
the plate begins to slide significantly. The layer begins to dilate while the applied stress continues to increase along 
BC, reaches a maximum at C, and then decreases along CD to its asymptotic value at E. The total variation of the 
vertical position of the plate is only a few micrometers (typically 5 //m) . 
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2. Layer dilation 



We find that the vertical position of the plate tends to its asymptotic value roughly exponentially with the sliding 
distance as shown in Fig. [?J and that the dilation h{x) is roughly independent of the driving velocity V (Fig. ||). This 
result suggests that the dilation rate dh/dt may be expressed as a function of the dilation h and of the plate velocity 
dx/dt as follows: 

dh h dx 

~dt = ~R dt 5 ( ' 

where R is the characteristic distance over which the layer dilates. The experimental behavior of dh/dt as a function 
of h and dx/dt is shown in Fig. |j. Although the law given in Eq. (||) is only roughly satisfied, the mean slope of the 
curve leads to R ~ 59 /im, which is approximately the mean radius of the glass beads. The corresponding exponential 
variation of h(x) shown in Fig. ^ demonstrates good agreement between Eq. (||) and the experimental data. 

We measure no systematic dependence of the total dilation Ah = ft(oo) — h(0) on the driving velocity V for layers 
prepared in the same way. Any variation is, at most, 1 (im over the full velocity range accessible to the experimental 
setup. Since we expect the initial compaction of the layer to be the reproducible (given identical preparation), 
we conclude that the mean density of the sheared granular layer does not depend on the shear rate within the 
experimental resolution. In contrast, the total dilation Ah decreases when the normal applied stress is increased; Ah 
typically decreases by 0.3 /im when the mass M of the plate is increased by 1 g. 

The total dilation Ah depends strongly on the initial conditions. For instance, we show in Fig. [l(] the different 
behaviors of h and d as functions of time in two cases: (1) The horizontal stress is released prior to the experiment. 
[The layer is prepared as described in Sec. III. A with results shown in Fig. ||.] In this case we find Ah ~ 5 /im. 
(2) The horizontal stress is not released between runs. [The plate is initially pushed at constant velocity (20 /im/s) 
until the steady-state regime is reached. The motion of the translator is suddenly stopped and the plate stops at a 
well-defined horizontal applied stress (F = 3.2 x 10 -2 N ~ Fd). The translator motion is then started again.] In this 
second case we find that the total dilation of the layer during the motion is only Ah ~ 1 /im. The smaller dilation 
observed in case (2) suggests that the continuously applied horizontal stress prevents the layer from compacting freely 
between runs. 

The small magnitude of the total dilation, roughly 10% of the bead radius, implies that the shear zone is localized 
to at most a few layers of beads. The near constancy of the vertical displacement during sliding suggests that shear 
may organize the beads into horizontal layers. 

In summary, a significant dilation of the granular layer accompanies the horizontal motion of the plate. During the 
transient regime, the layer dilates from its initial compaction to allow the horizontal motion of the plate. This slight 
expansion of less than one tenth of the bead diameter occurs over a characteristic horizontal displacement comparable 
to the bead radius. The total dynamical dilation decreases when the normal stress is increased and does not depend 
on the driving velocity. Experiments performed with 200 and 500 /im diameter beads show that the total dilation 
scales with the bead size. 



3. Frictional force 

The frictional force F reaches a maximum F max during the initial transient while the layer dilates (label C in 
Fig. ^). For a given layer under the same experimental conditions, the measured value of F max is reproducible to 
within 10%. As explained in the following, this experimental scatter originates essentially in fluctuations of the initial 
compaction. This scatter is small enough not to interfere with a measurement of the maximum frictional force as a 
function of V and Ah. 

In contrast to the behavior of the steady-state frictional force Fd, the maximum frictional force F max increases 
with V over the whole range of accessible driving velocities. Nevertheless, we find that F max increases only slowly for 
< 100 /im/s as shown in Fig. [□]. 

When layers of new material are prepared, the initial compaction fluctuates somewhat. We infer this fact from 
the dependence of the maximum frictional force F max on the total dilation Ah for many layers that have not been 
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subjected to an initial horizontal stress. The results are plotted in Fig. |l2], where one can see that the maximum 
frictional force increases linearly with Ah: 

F m ax = F d + a + (3 Ah, (4) 

whith a = (7.9 ± 0.4) 1(T 3 N and j3 = (8.4 ± 1.0) N//xm. 

Thus, the initial overshoot of the frictional force is at least partially related to the additional energy the system 
requires to dilate. In Fig. [l3], we show the variation of the frictional force F with the dilation rate dh/dt during a 
single transient event. As soon as the plate moves significantly in the horizontal plane (B and afterwards), F increases 
roughly linearly with the dilation rate dh/dt and can be described by 

„ „ dh 

F^F d + v—, (5) 

with v = (10.7 ± 0.5) 10 3 Kg/s. Note that F depends linearly on the dilation rate dh/dt rather than on the dilation 
h itself. This result is surprising. If we suppose that the overshoot of the frictional force originates in the potential 
energy used to lift the plate, one would expect the frictional force to depend linearly on h: The energy balance between 
the additional potential energy acquired by the weight per unit time U = Mg (dh/dt) and the power provided by the 
driving system P = F (dx/dt) would lead to F = F d + Mg(dh/dt)/(dx/dt) = F d - Mgh/R according to Eq. (|). 
The experiment performed by Marone et al. at large stress agrees with this prediction jjj. Nevertheless, such a 
dependence is not observed in our experiments. The excess work done by the driving system during the overshoot is 
experimentally about twice the increase in the potential energy. 

We also notice from Fig. |l3| that Eq. (||) fails to describe the initial stage of the motion, when the velocity and the 
displacement of the plate remain small (typically dx/dt < 3 ^m/s and x < 3 fim). The total displacement of the plate 
is then only 6% of the bead radius and the velocity is so small that we can expect the system to respond as an elastic 
medium. However, the experimental setup does not currently allow us to study in detail the very early response of 
the granular layer to a stress loading and unloading because of play in the driving system. 

In summary, the experimental results exhibit the important role played by the dilation of the layer on the frictional 
force F, which depends roughly linearly on the dilation rate dh/dt when the plate experiences a significant motion 
in the horizontal plane. As a consequence, the maximum value reached by the frictional force during the transient 
for a given driving velocity V depends linearly on the total dilation Ah of the layer between its initial state and the 
steady-state regime. 



D. Response to a static shear stress larger than the critical value 



Let us now consider the response of a compact granular layer to a static horizontal applied force F. The layer is 
prepared by driving the system in the steady-state regime at large velocity (typically 20 /im/s), stopping the motor 
suddenly, and pulling the spring back until the applied stress is fully released (F — 0). Afterwards, the spring is again 
pushed ahead but at lower velocity (typically 1 /zm/s) and is stopped at a given value of the applied stress F. 

We find that the granular layer can sustain the applied stress for a long time (several minutes) when F < 1.15 Fd 
typically. Nevertheless, the plate creeps slowly and the frictional force gradually declines. Small amplitude oscillations 
of the vertical position of the plate are observed (typically 0.2 /j,m in amplitude as shown in Fig. [l4|, case 1) without a 
mean rise of the plate; the horizontal velocity dx/dt (not shown) does not exhibit measurable oscillations. In contrast, 
when the applied stress is larger (e.g. F ~ 1.17 Fd, case 2 of Fig. |l4|), the granular layer sustains the stress only 
for a few minutes. The plate creeps slowly at first while a significant dilation of the layer is observed (typically 
1 /im). After a few minutes a large slip event occurs during which the stress is released. The plate slides rapidly for 
a few tens of micrometers in the horizontal plane and compacts by about 2 /im vertically. After the slip event, the 
horizontal motion of the plate is again hardly noticeable and vertical oscillations are observed as for lower applied 
stress. The plate can then sustain the remaining stress, which is now smaller, for a long time. We do not understand 
the oscillations of the vertical position of the plate, but they cannot be experimental artifacts, to the best of our 
knowledge. 
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IV. DISCUSSION AND CONCLUSION 



A. Recapitulation of the main experimental results and empirical model 



The experiments performed on immersed granular material allow us to propose an empirical model for the mechan- 
ical behavior of a sheared granular layer under very low stress. The model may also be applicable to dry materials 
but a stiffer measurement system would be required to verify this. The main results are as follows: 



1. Layer dilation 



Any horizontal motion of the plate involves a dilation of the granular layer: 

(i) The vertical position of the plate h tends roughly exponentially (Fig. to its asymptotic value over a distance 
R, which is approximately the bead radius. The dilation rate dh/dt then obeys Eq. (|3|). 

(ii) The total dilation Ah of the layer does not depend on the driving velocity V when the plate slides continuously. 

(iii) The total dilation Ah of the layer decreases when the normal applied stress is increased. 

(iv) The total dilation Aft, of the layer scales like the bead size. 



2. Frictional force 



After a transient, the plate generally slides continuously. The measured frictional force in the steady-state regime 
is proportional to the normal applied stress (Fig. ||) and does not depend on the driving velocity V over the whole 
velocity range accessible to the experimental setup (Fig. ||). Nevertheless, the existence of a stick-slip motion (Fig. ||) 
and the response of the granular layer to a static applied stress (Fig. |l4|) are consistent with a velocity weakening of 
the granular layer at very small velocity (typically 0.1 /jm/s). 

During the transient, while the layer is dilating significantly, the frictional force F depends roughly linearly on the 
dilation rate dh/dt (Eq. (jo] ) and Fig. and the frictional force reaches a maximum value that increases with the 
driving velocity V (Fig. and increases linearly with the total dilation Ah of the layer (Eq. (0) and Fig. O). 



3. Empirical model 



The motion of the plate is approximately governed by the differential equation given by Eq. (Q) in which F must 
be replaced by its empirical expression proposed in Eq. |^. The vertical position of the plate obeys Eq. ^. The system 
of differential equations that governs the time evolution of x and h can be written: 

, r d 2 x , , T _ . „ dh 
M w = k(Vt- X )-F d -v- 

dh h dx 

~dl ~ ~R~dt ^ ' 



with the initial conditions: 



x(0) = ; §(0) = 

MO) = -ho. (7) 



The initial value — ho of h is provided by the experiment. The system of differential equations [Eqs. (g) and ([?])] is 
then integrated numerically using the Runge-Kutta method jfj| . 

In the next section, we discuss the results of our simplified empirical model and compare them to the experimental 
measurements. 
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4- Comparison of the empirical model to the experimental results 



We show in Fig. 15 the result of integrating the empirical model and comparing it to the experimental data. The 
model is expected to describe the dynamics of the plate only from B to E. We take as the initial value of h its value 
at B. The model correctly describes the dilation of the layer. However, the agreement between the experimental and 
the theoretical instantaneous values of the force is imperfect; the small discrepancy is mainly due to the fact that the 
model Equation (||) does not account for the detailed behavior of dh/dt in the cycle shown in Fig. @. Nevertheless, 
the empirical model is in qualitative agreement with the experimental results and allows one to predict an increase of 
the maximum frictional force F max with the total dilation Ah and the driving velocity V. 

On the other hand, the empirical model cannot describe the very early stage A to B of the plate motion and the 
response of the granular layer to a static applied stress. Further experimental studies would be required to obtain a 
complete description of these phenomena. Our experimental setup does not allow us to perform these experiments in 
its present configuration. 



B. Comparison to previous work 



Both this work and previous studies []3| point out the important role played by the dilation on the frictional 
properties of the sheared granular material. However, the results obtained at small stress differ significantly from 
previous experimental results obtained at much larger stress. 

To explain these differences, we first summarize the friction law given by Marone et al. H. The first law decribes 
the dependence of the frictional force on the slip velocity V and on the surface's slip history via a state variable \t 

F = Mg Uo + + a In , (8) 

where F is the shear stress, Mg is the normal applied stress, fj,Q is a constant (which can be understood as the overall 
frictional coefficient), V* is an arbitrary reference velocity, and a and b are two empirical constants. The evolution of 
the state variable is governed by 

— = * + In — , (9 

dt R\ V* ) y ' 

where R is a characteristic distance over which the frictional force changes following a change in the slip velocity V. 
According to Eqs. (||) and (J)|), the steady-state value of the frictional force F d depends on V and dF d /d(\n V) — a — b. 
The experiments performed at large normal applied stress qualitatively agree with this theoretical description. 

We find experimentally that the steady-state value of the frictional force Fd does not depend on the driving velocity 
V at low stress. In the absence of any velocity strengthening or weakening, a = b so that dFd/d(ln V) = 0, and the 
frictional force reads 

F = F d - aMgR—^, (10) 

where we have set Mgfi = Fd because d^/dt = in the steady-state regime. Eq. |l(] agrees with the energetic 
argument (equality of work done and potential energy gain) mentioned in Sec. III. C. 3 provided that = —h/R. 
We tentatively assume this connection between the state variable \& and the vertical displacement h in our work in 
order to allow an interpretation of our results in term of their theory. With this assumption the model of Marone et 
al. leads to F = Fd + Mg^. Our experiments performed at small applied stress disagree qualitatively with these 
results. Indeed, the frictional coefficient is found to increase linearly with the dilation rate dh/dt (Eq. |5|) rather than 
with dh/dx. However, in making a comparison one should also note that the normal and tangential forces are applied 
independently to the granular material in our work, while in Ref. , the two are equal. 
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C. Conclusion 

The immersion of the material presents several experimental advantages. First, it allows one to work in well-defined 
conditions and to suppress any variability related to humidity changes. Efforts to eliminate water are rarely adequate 
because of adsorption. Second, a continuous motion of the sliding plate is observed even at low driving velocity (that 
is, stick-slip motion is avoided) and an extremely precise study of the steady frictional properties of the granular layer 
is then possible. The transient behavior that precedes the steady-state continuous motion allows a precise study of 
the dynamics of the frictional force for granular materials under low stress. A comparison with previous work shows 
that there are significant differences from the high stress case most relevant to geophysics. 

Certainly, there are different physical processes at work at high pressures, where the individual particles can be 
fractured by the stress and plastic flow may also occur. We hope to obtain additional insight by imaging the granular 
layer during the motion of the plate. The response of the granular layer on very long time scales to static applied stress, 
including the slow strengthening of the material in the presence of stress, will be analyzed in a further publication [fl4|| . 
We believe that our proposed friction law would also apply to dry material at low stress, although further experiments 
would be required to demonstrate this. 
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FIG. 1. Schematic diagram of the experimental setup for studying a fluid-saturated granular layer with sensitive measurement 
of horizontal and vertical positions. 



FIG. 2. Typical behavior (a) of the spring displacement d(t) and (b) of the vertical position h(t) as a function of time t 
(k = 189.5 N/m, M = 14.5 g, V = 28.17 /jm/s). [The different stages of the plate motion, labelled A-E, are discussed in section 
III.C.l.] 

FIG. 3. Position of the plate a; as a function of time t (k = 189.5 N/m, M = 14.5 g, V = 28.17 /xm/s). The dashed line 
indicates the asymptotic regime. 
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FIG. 4. Frictional force F in the steady state dynamical regime as a function of the mass M: The slope gives fid = 0.236±0.004 
(k = 189.5 N/m, V = 28.17 pm/s). The circle denotes the plate used in most of the experiments. 

FIG. 5. Frictional force F in the steady-state regime as a function of the pulling velocity V. There is no evidence of a 
dependance of the frictional coefficient fid on V (k = 189.5 N/m, M = 14.5 g). 

FIG. 6. Typical stick-slip motion observed for very small velocities. The rising parts of d(t) correspond to dx/dt = 0; The 
plate is at rest and the applied stress increases linearly with time. The sudden decreases of d correspond to the slip events; The 
maximal velocity is then about 2 /im/s. The amplitude of oscillation is about 10 /im around the mean value < d >= 192 /im 
(k = 189.5 N/m, M = 14.5 g, V = 0.11 /im/s). 

FIG. 7. Vertical displacement ft as a function of the horizontal displacement x of the plate (k = 189.5 N/m, M — 14.5 g, 
V = 28.17 (im/s). Dots - experimental points. Line - exponential interpolation with R = 59 fim in Eq. ^. The slope is finite 
at x = 0. 

FIG. 8. Vertical position h(x) as a function of the horizontal position x of the plate for different velocities. The granular 
layer dilates over a distance comparable to the bead radius (k — 189.5 N/m, M — 14.5 g). 



FIG. 9. Dilation rate dh/dt as a function of h{dx/dt). The linear interpolation of Eq. g leads to R ~ 59 fim, which is 
comparable to the bead radius (k = 189.5 N/m, M = 14.5 g, V = 28.17 /im/s). The arrows indicate increasing time. 



FIG. 10. Behavior (a) of the spring displacement d(t) and (b) of the vertical position h(t) as functions of time t in two 
different cases: (1) The horizontal stress is released before the experiment; (2) The horizontal stress is continuously applied. In 
the second case, the layer is initially less packed and, as a consequence, the total dilation Ah observed during the experiment 
is less (k = 189.5 N/m, M = 14.5 g, V = 28.17 /im/s). 



FIG. 11. Maximum frictional force F max as a function of the driving velocity V (k = 189.5 N/m, M = 14.5 g). 



FIG. 12. Maximum frictional force F max as a function of the total dilation Ah = h(oo) — h(0). The layer is initially unstressed 
(k = 189.5 N/m, M = 14.5 g, V = 28.17 /im/s). 

FIG. 13. Frictional force F as a function of the dilation rate dh/dt showing that F increases roughly linearly with dh/dt 
between B and D (k = 189.5 N/m, M = 14.5 g, V = 28.17 /tm/s). The initial oscillations from A to B are experimental 
artifacts due to the differentiation of experimental data containing noise. 

FIG. 14. (a) Applied stress F and (b) vertical position h as functions of time when a static stress larger than the critical 
value is applied (k = 189.5 N/m, M = 14.5 g, V = 0). Two cases are shown: (1) F ~ 1.1 F d ; (2) F ~ 1.17 F d . 

FIG. 15. Experimental behavior compared with the theoretical model (a) of the spring displacement d(t) and (b) of the 
vertical position h(t) as a function of time t (k = 189.5 N/m, M = 14.5 g, V = 28.17 /im/s). dots - experimental data, lines - 
empirical model. 
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FIG. 1. Schematic diagram of the experimental setup for studying a fluid-saturated granular layer with sensitive measurement 

of horizontal and vertical positions. 
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Typical behavior (a) of the spring displacement d(t) and (b) of the vertical position h{t) as a function of time t 
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FIG. 5. Frictional force F in the steady-state regime as a function of the pulling velocity V. There is no evidence of a 
dependance of the frictional coefficient /id on V (k = 189.5 N/m, M = 14.5 g). 
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FIG. 6. Typical stick-slip motion observed for very small velocities. The rising parts of d(t) correspond to dx/dt = 0; The 
plate is at rest and the applied stress increases linearly with time. The sudden decreases of d correspond to the slip events; The 
maximal velocity is then about 2 /im/s. The amplitude of oscillation is about 10 /im around the mean value < d >— 192 /im 
(k = 189.5 N/m, M = 14.5 g, V = 0.11 (j,m/s). 
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FIG. 7. Vertical displacement ft as a function of the horizontal displacement x of the plate (k — 189.5 N/m, M = 14.5 g, 
V — 28.17 fim/s). Dots - experimental points. Line - exponential interpolation with R = 59 /im in Eq. ^. The slope is finite 
at x = 0. 
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FIG. 8. Vertical position h(x) as a function of the horizontal position x of the plate for different velocities. The granular 
layer dilates over a distance comparable to the bead radius (k — 189.5 N/m, M = 14.5 g). 
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FIG. 9. Dilation rate dh/dt as a function of h(dx/dt). The linear interpolation of Eq. ^ leads to R ~ 59 /im, which is 
comparable to the bead radius (k = 189.5 N/m, M = 14.5 g, V = 28.17 /im/s). The arrows indicate increasing time. 




FIG. 10. Behavior (a) of the spring displacement d(t) and (b) of the vertical position h(t) as functions of time t in two 
different cases: (1) The horizontal stress is released before the experiment; (2) The horizontal stress is continuously applied. In 
the second case, the layer is initially less packed and, as a consequence, the total dilation Ah observed during the experiment 
is less (k = 189.5 N/m, M = 14.5 g, V = 28.17 /jm/s). 
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FIG. 11. Maximum frictional force F max as a function of the driving velocity V (k = 189.5 N/m, M = 14.5 g). 
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FIG. 12. Maximum frictional force F ma x as a function of the total dilation Ah — h(oo) — h(0) . The layer is initially unstressed 
(k = 189.5 N/m, M = 14.5 g, V = 28.17 //m/s). 
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FIG. 13. Frictional force F as a function of the dilation rate dh/dt showing that F increases roughly linearly with dh/dt 
between B and D (k = 189.5 N/m, M — 14.5 g, V = 28.17 /im/s). The initial oscillations from A to B are experimental 
artifacts due to the differentiation of experimental data containing noise. 
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FIG. 14. (a) Applied stress F and (b) vertical position h as functions of time when a static stress larger than the critical 
value is applied (k = 189.5 N/m, M = 14.5 g, V = 0). Two cases are shown: (1) F ~ 1.1 F d ; (2) F ~ 1.17 F d . 
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FIG. 15. Experimental behavior compared with the theoretical model (a) of the spring displacement d(t) and (b) of the 
vertical position h(t) as a function of time t (k = 189.5 N/m, M = 14.5 g, V = 28.17 /im/s). dots - experimental data, lines - 
empirical model. 



